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Abstract 

It is known that for multi-level time-dependent quantum systems one can construct 
superadiabatic representations in which the couphng between separated levels is ex- 
ponentially small in the adiabatic limit. For a family of two-state systems with real- 
symmetric Hamihonian we construct such a superadiabatic representation and explicitly 
determine the asymptotic behavior of the exponentially small coupling term. First order 
perturbation theory in the superadiabatic representation then allows us to describe the 
time-development of exponentially small adiabatic transitions. The latter result rigor- 
ously confirms the predictions of Sir Michael Berry for our family of Hamiltonians and 
slightly generalizes a recent mathematical result of George Hagedorn and Alain Joye. 

1 Introduction and main result 

The decoupUng of slow and fast degrees of freedom in the adiabatic limit is at the basis 
of many important approximations in physics, as, e.g., the Born-Oppenheimer approxima- 
tion in molecular dynamics and the Peierls substitution in solid state physics. We refer to 
[BMKNZI ITe| for recent reviews. Generically the decoupling is not exact and a coupling 
which is exponentially small in the adiabatic parameter remains. However, this small cou- 
pling has important physical consequences, as it makes possible, e.g., non-radiative decay 
to the ground state in molecules. Since Kato's proof from 1950 Ka the adiabatic limit of 
quantum mechanics was considered also as a mathematical problem, with increased activity 
during the last 20 years. Some of the landmarks are |Nei| lASYl | JoPfi] |Ne2[ IBaJoj . 

We consider a two-state time-dependent quantum system described by the Schrodinger 
equation 

{iedt - H{t))^{t) = (1) 
in the adiabatic limit £ — > 0. For the moment we take the Hamiltonian H{t) to be the 



real-symmetric 2 x 2-matrix 



H{t) - p{t) 



cosd{t) sm0{t) 
sme{t) -cose{t) 



The eigenvalues of H{t) are ±p(i) and we assume that the gap between them does not close, 
i.e. that 2p(t) > .g > for aU t e M. 

As to be explained, even for this simple but prototypic problem there are open mathemat- 
ical questions. In order to explain the concern of our work, namely the time-development of 
the exponentially small adiabatic transitions, let us briefly recall some important facts about 
QJ. Let Uo{t) be the orthogonal matrix that diagonalizes H{t), i.e. 

_ / cos(0(i)/2) sin(0(i)/2) \ 
^°^^> ~ y smi9{t)/2) ^cos{e{t)/2) J ■ 

Then the Schrodinger equation in the adiabatic representation becomes 

Uo{t) {iedt - H{t)) Uo{t)m {iedt - H^{t))^-{t) = 

with 

/ Pit) i§9'{t) \ 

^■•'" = (-W(o %«) j *''" = ™- 

Here and henceforth, primes denote time derivatives. First order perturbation theory in 
the adiabatic representation (cf. proof of Corollary ^ and integration by parts yields the 
adiabatic theorem |BoFo[ IKa| : The off-diagonal elements of the unitary propagator K^(t, s) 
in the adiabatic basis, i.e. the solution of 

iedtK^it, s) - HI{t)K^{t, s) , K^{s, s) = id , 

vanish in the limit e ^ 0. More precisely, let 

--(Jo). ---ill)' 

which project onto the adiabatic subspaces in the adiabatic representation. Then 

\\P^K-{t,s)P+\\^Oie). (4) 

Therefore the transitions between the adiabatic subspaces are 0(e). This bound is optimal in 
the sense that in regions where 9(t) is not constant the leading order term in the asymptotic 
expansion of P_ K^{t, s) P+ in powers of e is proportional to e. 

However, if \mit^±oc 9' {t) = then in the scattering limit the transitions between the 
adiabatic subspaces are much smaller: if the derivatives of 6* e C°°(M) decay sufhciently fast, 
then for any n G N 

A{e) := hm || F_ K^,[t, -t) P+\\ = . (5) 

r — >oo 

If 9 is analytic in a suitable neighborhood of the real axis, then transition amplitudes are 
even exponentially small, .4(e) — 0{e~'^/'^) for some constant c depending on the width of 
the strip of analyticity, see |JoPfi| IMaj . 
It is well understood, see ILellGal lNei 



how to reconcile the apparent contrariety between 
the smallness of the final amplitudes in ((SJ and the optimality of Q: the adiabatic basis is 
not the optimal basis for monitoring the transition process. For any n g N there exist unitary 
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transformations U"{t) such that the Hamiltonian in this n}^ superadiabatic representation 
takes the form 

H^'it) = ( ^li^l ) with p-{t) = Pit) + 0{s') and \c-{t)\ = 0{s-+') . (6) 



In the n^^ superadiabatic basis the off-diagonal components of the propagator and hence also 
the transitions are of order ©(e"), i.e. there are constants C„ such that 

\\P^K-{t,s)P+\\<CnS-. (7) 

In the scattering regime, where d(t) becomes constant, the superadiabatic bases agree with 
the adiabatic basis, i.e. \imt-,±oo U" (t) = Uo{t), and therefore the bound in (Q basically 
yields ll^Jl. Typically lim^^oo C„e" = oo for all e > 0, i.e. choosing n larger while keeping e 
fixed does not necessarily decrease the bound in Q. However, one can choose — n{e) in 
such a way that C„^£"^ is minimal. If 9 is analytic, one obtains the improved estimate 

II P-K^'it,s)P+\\ = ©(e-^/'^) 

in the optimal superadiabatic basis n^, see |Ne2l poPf2| . 

More interesting than bounds on .4(e) is its actual value. Since A{e) is asymptotically 
smaller than any power of e, this question is beyond standard perturbation theory. For the 
case of analytic coupling 0, asymptotic formulas of the type 

yt(e) = Ce-^ (l + 0(e)) (8) 

have been established, see e.g. |JKPI Uo) . where the constants C and tc depend on the type 
and location of the complex singularities of 9'{t)/p{t). However, these results are obtained 
by solving |^ not along the real axis but along a Stoke's line in the complex plane. As 
a consequence they give no information at all about the way in which the exponentially 
small final transition amplitude A{s) is build up in real time. This question of adiabatic 
transition histories is the concern of our paper. Berry Be^ and, in a refined way Berry and 
Lim I BeLiI ILiBe| , gave an answer on a non-rigorous level and explicitly left a mathematically 
rigorous treatment as an interesting open problem. Only very recently Hagedorn and Joye 
|Ha,Toj succeeded and confirmed Berry's results rigorously for a specific Hamiltonian. 

Although our work has been strongly motivated by the findings of Berry, our approach 
is slightly different. Let us first state our main result before we discuss its relation to the 
earlier ones. Without loss of generality we assume that p{t) = ^- It was observed in |Be| . 
that this can always be achieved by transforming |^ to the natural time scale 



r(i) ^2 f\ 
Jo 



g{s) ds . 

However, we can only treat a rather special class of Hamiltonians, since we must assume that 
in the natural time scale the coupling has the form 



(9) 



with 7 G M and tc > 0. In other words we assume 

^W^K"n'S »W^^7«c.a„(i|^ ,10) 
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We shall comment below on the meaning of this special choice and remark here that the 
Hamiltonian in jHa.Toj is with 7=5- 

Our main result is the construction of an optimal superadiabatic basis in which the 
coupling term in the Hamiltonian is exponentially small and can be computed explicitly at 
leading order. This optimal basis is given as the n**^ superadiabatic basis where < cr^ < 2 
is such that ^ 

rig = — — 1 + (jg is an even integer. (11) 

£ 

Theorem 1. Let H{t) he as in and as in and let £q > be sufficiently small. 

Then for every e € (0,£o] one can construct a family of unitary matrices U"'^{t) G C^^^, 
depending smoothly on t € such that 

\\U-^it)-Uom=o(^j^^ (12) 

and 

U-^it) (isd, - Hit)) U-^*it) = iedt ^l^^ ) . (13) 



H-^it) 



Here 



and for every a < | 



P^^{t) = l + o^ ' 



2 \l + t^ 



2 



2e . /7r7\ _tc_ ft i . cTst 



cTH) -2i^- sin(-j e-- e"-. cos ^ - - ^ + ^ ) +O(0"(M)) , (14) 
with 

( e-expf-^fl + i^)) tf\t\<t,, 
(l>"{e,t) = \ ' (15) 

Remark 1. The explicit term in is of order ©(c^*"/^) only for times |t| = 0{^/e). For 
larger times all terms in c"^ are exponentially small compared to the leading exponential 
Q-tc/e ^ As a consequence, Taylor expansion of the cosine in c"^ around t/e for \t\ = 0(-\/e) 
shows that it can be replaced by cos(i/e) at the cost of lowering a to a < 1: for every a < 1 

2e . /7r7\ _tc / 1 
sin I e e 



c^' (t) = 2i \/ — sin j e-- e"^ cos j + O (0" (e, t)) . 

Remark 2. The slow time decay of the error in (|15|l for large times is due to the fact that 
rig is optimal for t near 0, but not for large t. 

Remark 3. Taking defined in pi|l odd instead of even would yield slightly different 
off-diagonal elements in the effective Hamiltonian H^^(t). However, the resulting unitary 
propagator, cf. Corollary^ would be the same at leading order. See the end of Section 5 for 
a discussion of this somewhat surprising fact. 

Let us shortly explain the idea of the proof of Theorem ^ and at the same time the 
structure of our paper. First we construct the n^^ order superadiabatic basis as in ((BJ in two 
steps: in Section 2 we construct the projectors on the superadiabatic basis vectors and in 
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Section 3 we construct the unitary basis transformation J7" (t) . We cannot use existing results 
here, e.g. j(Tal|Ne2l , since we need to keep carefuUy track of the exact form off the off-diagonal 
terms c" (i) of the superadiabatic Hamiltonian, and since we aim at a scalar recurrence relation 
instead of a matrix recurrence relation for the c"(i)'s. The main mathematical challenge is 
the asymptotic analysis of the resulting recurrence relation, which is done in Section 4. This 
is also the only part where we have to assume the special form Q for 6'. Theorem ^ then 
follows by choosing the order n of the superadiabatic basis as in , a choice which minimizes 
c"(i) near t = 0. The details of this optimal truncation procedure and the proper proof of 
Theorem^are given in Section 5. Finally in Section 6 we use first order perturbation theory 
in the optimal superadiabatic basis in order to obtain the following Corollary, in which we 
abbreviate 

A{t, s) := arctan(t) — arctan(s) . 

Also recall that erf: M (—1, 1) with erf(a;) = lo '^^^ da; switches smoothly and mono- 
tonically from erf(— oo) = —1 to erf(oo) — 1. 



Corollary 1. The unitary propagator in the optimal superadiabatic basis 

K-'it,s) = 

i.e. the solution of 



k+{t,.3) k^{t,s) 

k^{t,s) k-{t,s) 



iedtK^^ {t, s) = ff^' {t)K^' {t, s) , Kl'^ (s, s) = id , 

satisfies 



and 



/cf (i,s) = e=FT^ +0(eA(i,s)) (16) 



/ 7r7 \ _ to _ 



k^{t, s) = sin l^— ) e s e ( erf ( ) — erf 



0(V~ee-^Ait,s)) . (17) 



Outside the transition region, more precisely for \t\ > and \s\ > for some (3 < i, j j7| j 
holds with the error term replaced by 0{e"e'~'^ A{t, s)) for every a <1. 

Corollary n immediately implies the existence of solutions to |^ of the form 

( sin(^)e-¥el(erf(^)+l) ) + ^(^^"") ' ^''^ 

They start at large negative times in the positive energy adiabatic subspace and smoothly and 
monotonically develop the exponentially small component in the negative energy adiabatic 
subspace in a -^/e-neighborhood oi t = 0. Berry and Lim jBel IBeLi] argue that this behavior 
is universal: whenever 9' has the form 

9'{t) = ^^Y^ +0{\t±itc\"') forsomea>-l 

near its singularities ±itc closest to the real axis, then (|18f) should hold. For the Landau-Zener 
Hamiltonian, which describes the generic situation, one finds 7 = ^ and a = Hagedorn 
and Joye |HaJO| proved (|18|> for the Hamiltonian H1U|) with 7 = ^ • In the approach of Berry 
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and, slightly modified, of Hagedorn and Joyc, the optimal superadiabatic basis vectors are 
obtained through optimal truncation of an asymptotic expansion of the true solution of |^ 
in powers of e. 

In contrast, in our approach the optimal superadiabatic basis is constructed by approx- 
imately diagonalizing the Hamiltonian. The main advantage of "transforming the Hamil- 
tonian" over "expanding the solutions" is that the former approach can be applied, at 
least heuristically, to more general adiabatic problems, cf. |Te| . as for example the Born- 
Oppenheimer approximation. While we cannot control the asymptotics for the Born-Oppen- 
heimer model rigorously yet, the heuristic application of the idea yields new physical insight 
into adiabatic transition histories and new expressions for the exponentially small off-diagonal 
elements of the 5- matrix for simple Born-Oppenheimer type models, cf. |BeTe| . Therefore 
we see the rigorous results obtained in this paper also as a first attempt to justify the appli- 
cation of analogous ideas to more complicated but also more relevant systems. Furthermore, 
the concept of an adiabatically renormalized Hamiltonian was used to derive a criterion for 
selecting possible transition sequences in multi-level problems |WiMoj . 

For the specific problem the knowledge of two linearly independent solutions is of 
course equivalent to the knowledge of the propagator and the effective Hamiltonian in the 
optimal superadiabatic basis. Therefore we shortly explain which aspects of our result con- 
stitute an improvement compared to jHa.Toj : Most importantly, our proof does not rely on 
the a priori knowledge of the scattering amplitude A{e) . Indeed, our result yields for the first 
time a proof of © based on superadiabatic evolution, as expressed in Corollary|21 Moreover, 
we allow for a slightly larger class of Hamiltonians and obtain more detailed error estimates, 
which, in particular, give rise to close to optimal error bounds in the expansion of the S- 
matrix, cf. Corollary|21 Finally, we also get explicitly the next order correction in H17|l resp. 
H18(l . cf. Section 6. It should be noted, however, that the improved error estimates and the 
next order corrections could have been obtained also based on the proof in iHaJol . 

We finally turn to the scattering limit. Let K^{t, s) denote the propagator in the original 
basis and define the scattering matrix in the adiabatic basis by 

SI := hm Uo{t) K"^{t, -t) U*{-t) e^ , where Hq = ( ^ 



t- 



Since, according to (|12|) . for large negative and positive times the optimal superadiabatic basis 
agrees with the adiabatic basis, S'^ can be computed with help of the optimal superadiabatic 
propagator from Corollary 

Corollary 2. For (3 <1 we have 

( l + 0{e) 2sin(^) e-^(H-0(e^)) \ 

' 1^ 2 sin(^) e-^ (1-f 0(e/3)) l + 0{e) )' 

Proof. According to H12|) we have 

S"^ = lim e^ U^' (t) K°{t, -t) U','^ *{-t) e^ = lim e^ K"^^ {t, -t) e^ . 



Now the claim follows from inserting (|16l) and H17() with the improved error estimate outside 
of the transition region. □ 

From Corollary |2] we conclude that the transition amplitude is given by 

A{e) = \\P-SlP+\\ = 2 sin e"^ (l -I- ©(e^^)) , for any /3 < 1 , 
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which agrees with the results of [Jo], as explained in |BeLi| . 

We conclude the introduction with two recommendations for further reading: The numer- 
ical results of Berry and Lim jLiBe| beautifully illustrate the idea of optimal superadiabatic 
bases and universal adiabatic transition histories. The introduction of the paper of Hage- 
dorn and Joye |HaJo| gives a slightly different viewpoint on the problem and, in particular, 
a short discussion on how exponential asymptotics for the Schrodinger equation fit into 
the broader field of exponential asymptotics for ordinary differential equations. 

Acknowledgements: We are grateful to Alain Joye and George Hagedorn for many helpful 
discussions. 



2 Superadiabatic projections 

For the present and the following section we assume that H{t) has the form but with 
some arbitrary 9 G C°°(]R). The first aim is to construct time-dependent matrices tt'"^ g 

(7r("))2-7r(") =C'(£"+i), (19) 



iedt - H, 7r(") 



= C'(£"+0. (20) 



Here, [A, B] = AB — BA denotes the commutator two operators A and B. Likewise, we will 
later use [A, i?]^ = AB + BA to denote the anti-commutator of A and B. Equation (|19|l 
says that tt^^"^ is a projection up to errors of order £"+^, while (|20|l implies that 7r("^(i) is 
approximately equi variant, i.e. 

K°{t, s) 7r(") (s) = 7r("' {t) K°{t, s) -f ©(s") . 

Recall the K^{t,s) is the unitary propagator for Hence 7r(")(<) is an almost projector 
onto an almost equi variant subspace. 

We construct tt^"^ inductively starting from the Ansatz 

n 

TT^"^ =Y,^ke''. (21) 

By (|10f) . H has two eigenvalues ±1/2. Let ttq be the projection onto the eigenspace corre- 
sponding to -1-1/2, and 7r^°) = ttq according to H21|l . It is easily checked that H19|) and (|20|l 
are fulfilled for n = 0. In order to construct 7r„ for n > 0, let us write G„(t) for the term of 
order £"+^ in i.e. 

(^(n))2 _ ^(„) ^ £«+lG„+i + 0(£«+2) . (22) 

Obviously, 

n 

Gn+1 = T^jT^n+l-j- (23) 

i=i 

Proposition 1. Assume that tt*^"' given by \21\) fulfills \iy\) and \2U\) . Then a unique matrix 
TTn+i exists such that Tr'"^"'^' defined as in \21]) fulfills 1^1 Uj) and \2U\) . TTn+i is given by 

TTri+l = Gn+1 — — Gn+lT^O - i k^i, TTq] . (24) 

Furthermore tt^ is off-diagonal with respect to ttq, i.e. 

TToTT^TTo = (1 - 7roX(l - TTo) = 0, (25) 
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and Gn is diagonal with respect to ttq, i.e. 

7roG„+i(l — TTo) = (1 — 7ro)G„+i7ro 



0. 



(26) 



Remark 4. The fact that the superadiabatic projections are unique answers the question 
raised in \Be\ to which extend the superadiabatic basis constructed there is uniquely deter- 
mined. 

Remark 5. Our construction can be seen as a special case of the construction in jEmWe| . It 
was applied in the same context in [PSTllTte) . The role and the importance of the superadia- 
batic subspaces as opposed to the superadiabatic evolution have been emphasized by Nenciu 
|Ne2| . He constructs the superadiabatic projections for much more general time-dependent 
Hamiltonians. However, Nenciu's construction is less suitable for the explicit computations 
we need to perform. 



Proof. Let 7r("+i) be given by H21() and suppose tt*-") fulfills H19() and (|20|l . Let TTn-i-i be an 
arbitrary matrix, and define 7r("+i) = tt*^") + £"+^7f„+i. Then 

(^(n+l))2_^(„+l) ^ (^(n))2 



Using (|22|l . we see that terms of order vanish if and only if 

Gn+l — 7r„+l — [tTo, 7r„+l]^ = (1 ~ 7ro)7rn+l(l - TTq) — TroTi-n+lTTo. 



(27) 



Multiplying lf77|) with (1 — ttq) and with ttq on both sides and subtracting the results, we find 
that TTn+i must fulfill 



(1 — 7ro)7r„+l(l — TTq) + VTo^n+lTTo = Gn+l — [G„+i,7ro] 



(28) 



Similarly 



iedt- H,n^"+^'> = ieSf - i/, tt^"' 



e"+Mi£at-i?,*„ 



Again terms of order £"+^ vanish if and only if 



(29) 



Since ttq is the projector onto the eigenspace of H, we have ttqH = Httq = Ettq, where 
£^ = 1/2 is the positive eigenvalue of H, and similarly (1 — tto)H = H{1 — ttq) = ~E{1 — ttq). 
When we multiply H29|l first with with ttq from the left and with 1 — ttq from the right, then 
the other way round, and finally subtract the second result from the first, we get 



2i;(7ro7r„+i(l - ttq) + (1 - 7ro)7r„+i7ro) = -i [tt^, ttq] 
Now we divide ^ by 2E and add ^ to find 



G. 



n+l 



[G„+i, 7ro]_ 



2E 



(30) 



(31) 



Thus TTn+i is uniquely determined by the requirement that 7r("+i) should fulfill H19|) and H2U|I . 
On the other hand, [H, Gn+i — [G„+i, ttq]^] = and 

7I"0 K, TTo] TTo = (1 - TTo) [tT^, TTq] (1 - TTq) = , 

and thus 7r„+i given by the right hand side of H31I) indeed fulfills (|28() and H29|) . This shows 
existence. (j^HJ and now follow directly from l(?7|l and (Pi^ . □ 
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The calculation of tt^^"^ via the matrix recurrence relation H24I) and (|23() is now possible in 
principle, but extremely cumbersome. In order to make more explicit calculations possible, 
we introduce a special basis of R^^^. Recall that f7o(i) as defined in |5J) is the unitary 
transformation into the basis consisting of the eigenvectors of iJ, i.e. the adiabatic basis, 
and let ¥^{1) = g^U(){t). With P = P+ as in ^ we then have U§ = = id and 
PUqVqP = PVqUoP = 0, and ttq = UoPUq. Moreover, since Gi = by ^ imphes 

TTi = -^9'iVoPUo - UoPVo). (32) 

Motivated by this, we put 

X = VoPUo - UoPVo , Y = VoPVo - UoPUo , 

Z = VoPU„ + U„PVo , W = VoPVo + U„PUo . 

It is immediate that this is a basis of M^^^ for all t, and in fact 

Our reason for representing X through Z via Uq and Vq is that the following important 
relations now follow without effort: 

X' ==0, Y' = -e'z, z' = e'Y, (33) 

[X, Y]^ = [X, Z]^ = [Y, Z]^ = 0, -X^ = Y^ = Z^ ^W, (34) 
[X,no] = Z, [y,^o]-0, [Z,7To]=X, (35) 
W -[W,'Ko]^ = Y. (36) 

These relations show that this basis behaves extremely well under the operations involved in 
the recursion (|24|l . This enables us to obtain 

Proposition 2. For all n G N, 7r„ is of the form 

TTn = XnX + y„Y + Z„Z, (37) 

where the functions Xn , Un and z„ satisfy the differential equations 

x'„ = izn+i, (38) 

y'n = (39) 

= iXn+1 + O'yn. (40) 

Moreover, 

Xlit)^-'f'{t), y,{t)^z,{t)^0. (41) 

Remark 6. Hence, for all even n, a;„ = 0, while for all odd n, j/„ = z„ = 0. 

Proof. iflT)) was already noticed in . Now suppose 7r„ is given by ifTTI) . By (|^ and , 
Gn+i ~ [Gn+i, 7ro]_|_ is proportional to Y with a prefactor given through (HSl, and by ll^ . 
and 133l, 



7r„+i = '^{-XjXn+i^j + yjyn+i-j + ZjZn+i-j)Y + i{9'yn - z'^)X - ix'^Z. (42) 



Comparing with ||1T7|) shows and igUJl . To show , we use if^ . This gives 

Since ttoZttq — ttqXttq — and ttqYitq — ttq, the claim follows. □ 
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Remark 7. From H38|l through (|40() we may derive recm-sions for calculating Xn or z„, e.g. 

z„+2(i) = -^(4W+e'w(/ 9'it)z^it)dt + C^y (43) 
The constant of integration C must (and in some cases can) be determined by comparison 

with ggi- 

Using (|38|I -H40 |I . we can give very simple expressions for the quantities appearing in (|19|l 
and (|20|l . As for H20|l . we use H33|l and the differential equations to find 

k9t-i/,7r(")] = i£"+i< = i£«+i«X + (j/; + ^?'z„)y + (z;-0'y„)Z) = 

= -e"+\zn+lX +Xn+lZ). (44) 

Now we turn to (tt^"))^ — Tr'"^ the term by which tt'"^ fails to be a projector. Let us write 

n 

(^("))2 _ ^{n) ^ J2 e"+'=G„+i,fe. (45) 

k=l 

With our earlier convention, Gn+i,i — Gn+i- Explicitly, and (|45|l give 

71— k 

Gn+i,k = kfc, 7r„]_|_ + [TTfc+i, 7r„__i]^ + ■ ■ • ^ X/ (^6) 
Proposition 3. for each rt G N, t/iere ea;ist functions gn+i,ki k < n with 

((^("))2 _ ^ |^^^£«+^g„+, ,(t)j H^. (47) 

For each k < n, 

9n+l,k = 2i{xkZn+l - ZfeXn+i). 

Proof. By (I34II . each G„+i.fc is proportional to W. Using (|37|l additionally, we find [tt^, 7r,„]_|_ — 
2{-XkXrn + UkUrn + 2zfcZ„i)W^, and thus yields 

n — 

Qn+l^k — ^ ^ ^j+k'^n—j ~^ V j+kyn—j ~^ ^j+k^n—j- 
3=0 

Thus by using Proposition 12 

n — k 

g'n+l,k = X! K^j+k+lXn-j + Xj+kZn-j+l) - {0' Zj+kVn-j + O'y^+kZn-j) + 
j=0 

+0yj+kZn-j + O'Vj+kZn-j ~ K^j+k+lZn-j + Zj+k^n-j + l) = 
n — k 

= 2\{xkZn+l- ZkXn+l). 

The last equality follows because the sum is a telescopic sum. □ 

Since 14^ = id is independent of t. Proposition O gives the derivative of the correction 
(7r("))2 _ 7r(") to a projector. As above, this gives an easy way for estimating the correction 
itself provided we have some clue how to choose the constant of integration. 
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3 Construction of the unitary 

We now proceed to construct the unitary transformation U" into the n^^ superadiabatic 
basis. By (EU and lj23, tt^") is self-adjoint. Thus it has two orthonormal eigenvectors Vn 
and Wn- Let 

cos(6i/2) \ _ / sin(6l/2) 



sm{e/2) J' V -cos(0/2) 

be the eigenvectors of ttq , and write 

Vn — oiVo + f3wQ, Wn = awo — l3vo (a,/3eC). (48) 

We make this representation unique by requiring < a G M. Let [/" be the unitary operator 
taking w„) to the standard basis (61,62) of , i.e. 

[/;' = ei< + e2<, (49) 

where all vectors are column vectors. Note that the definition (O of Uq is consistent with 
for n = 0. U" diagonalizes tt^"^ , thus 



U A2 



(50) 



where Ai^2 are the eigenvalues of Tr*^"-*. Although a,f3 and Ai_2 depend on n, e and t, we 
suppress this from the notation. 

Lemma 1. 

Proof. The calculations are straightforward and we only show the second equality. First note 
that Wq = —Wo and = wq- Thus 

Ur' = ((a' + P)vo + {13' - a)wo)el + {{a - ^Vo + («' + ;^)^«o)6^, 
and using the orthogonality of uq and wq yields the claim, 

UoU^* — ei{a' + (3)el + ei(a — /3 )62 + 62(/3' — a)el + 62(0:' + (3)el . 

□ 

It will turn out that /3, a'a, and /3' are small quantities, A'i,A2, and A2 are even much 
smaller, while and Ai are large, i.e. of order 1 + 0{e). This motivates the form in which 
we present the following result. 

Proposition 4. Suppose Ai 7^ A2. Then for each n G 

/ 1 Q^e^+i _ \ \ 

[Xn+l — Zn+1 ) 



C/;(ie5,-iI)C/;*=iea,- I .^„,, 2 a,-a. v --^^ + i?, 

V A, -A, - Zn+l) "n / 



with 

R = 



eIm(/3(2a + /3')) + |/3p -xr^(^"+l + ^"+1) 
- -eIm(F(2Q + /3')) - 
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Proof. Let us write {iedt - H)U^ * = (M^j), i,j £ {1,2}. and M2^2 are calculated in 

a straightforward manner, using Lemmas together with the fact UqHUq ~ ^ -1/2 )• 



U^i^iedt - H)Ul'* = iedt + ieU^U^UoU'" * - U^U^UqHU^UqUI 

(3 a 



Carrying out the matrix multiphcation yields 



Afi^i = -Af2,2 = ie9t + ie{{a{a' + 13) + PiP' - a)) - ^{a^ - (51) 

We now use + ^ 1 to obtain = 2aa' + (3']3 + p' f3 = 2Re{aa' +]3P') and - |/3|2 = 
1 — 2|/3p. Plugging these into (|51|l gives the diagonal coefficients of M. Although we could get 
expressions for the off-diagonal coefficients by the same method, these would not be useful 
later on. Instead we use (|^ . i.e. U^*D — n^^^^U"* together with and obtain 

Ul\iedt - H)Ul'*D = DU^iiedt - H)U^* - e"+iC/;(z„+iX + Xr,+iZ)Ul'* . (52) 

By multiplying l(5^ with e^e* from the left and by e^e^ from the right (j, fc G {1,2}) and 
rearranging, we obtain 



(Afe - A, ) e, e* {iedt - H) *efc = (53) 
= -£"+'ej e* U'2{zn+iX + x,,+^Z) *ek el - i<5fc,, e X'^ e, e* . 



From the equalities UoXUq = ^ ) ' ^"^^0 ^ '^1 0^ ) ^'^^ Lemma^Jwe obtain 

The expressions for Mi, 2 and M2.1 follow by taking fc ^ j in (|53(l . □ 

We now use our results from the previous section to express a, (3 and A1.2 in terms of 
a;fe, Uk and Zk, k < n. Let us define 

e = e(",e,i) =ELie'^fcW> (54) 
V = V{n,e,t) ^Yl=i£^yk{t), (55) 
C = C(",e,i) =ELi£'^fcW- (56) 

Moreover, let 

g^.9(n,£,t) = ELie"+V.+i.fc(0 (57) 
be the quantity appearing in H47() . 

Lemma 2. T/ie eigenvalues of n^"^"^ solve the quadratic equation 

^1,2 - ^1,2 -5 = 0. 
Proof. By H5()|l and Proposition |2| we obtain 



^1 \ [/"((7r("))2 - ^(")) C/"* = U^gWU^* =( ^ ° 
\^UA2-A2y ■ \ Q 9 



□ 
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Lemma 3. 

a^(Ai - A2) = 1 - - A2, and a/?(Ai - A2) = - C- 
Proof. We have 

TT^"' = Aii;„< + A2W„w;. (58) 

Plugging in (|48(l . we obtain 

7r''"'wo = Aiaw„ - A2/3w„ = (Aia^ + A2|/3|^)wo + (Ai - A2)q;/3wo = 
= (a^(Ai - A2) + A2)wo + (Ai - A2)a/3wo- 

In the last step, we used |/3p + = 1. On the other hand, from and we have 

n 

TT^'''> =TTo + Y.^''('^kX + ykY + ZkZ), (59) 

k=l 

and since Xuq = Zvq — —wq, ttqVo ~ vq and yt;o — — Wo, we find 

^(")«o = (1 - r?)i;o - (e + C)wo. 
Comparing coefficients finishes the proof. □ 

Theorem 2. Let Eq > be sufficiently small. For e € (0,eo] assume there is a bounded 
function q onM. such that (,{t), r]{t), C(t) and their derivatives f'(i), ?7'(0i C'(^) "'^fi '^^^ bounded 
in norm by eq{t). Then 



U^iiedt- H)U^* ^ (60) 
/ i + 0(62^2) £"+i(a;„+i-z„+i)(l + 0(£g)) \ 

e"+H-^n+i - ^n+i) (1 + 0(eg)) -i + 0(e2g2) )• 



ie5f 



Proof. From H59(l and our assumptions it follows that 7r(")-7ro = ^^(eg). Thus A 1 = l + 0{eq) 
and A2 — 0{eq), and from LemmaElwe infer g = 0(e) and 

Ai = i (1 + VTT4^) , A2 = i (1 - 71+4^) . 

Since Ai — A2 7^ 0, Lemma |21 yields 

2 l + yrT4^-277 ^ -c-c 



Hence = 1 + 0(£g), and /3,/3' and aa' = (a^)'/2 are all 0(6(7). Plugging these into the 
matrix R in Proposition 0] shows the claim. □ 

4 Solving the recursion: a pair of simple poles 

In order to make further progress, we need to understand the asymptotic behavior of the 
off-diagonal elements of the effective Hamiltonian in the n*^ superadiabatic basis for large n. 
According to Ht)U|l this amounts to the asymptotics of a;„ and z„ as given by the recursion 
from Proposition 121 It is clear that the function 6' alone determines the behavior of this 
recursion. We will study here the special case 

0'(t) = - = .g^. 



13 



The reason lies in the intuition that the poles of 6' closest to the real axis determine the 
superadiabatic transitions, and that these transitions are of universal form whenever these 
poles are of order one, see |Bel IBeLi] for details. As in |Ha.To| . we have to restrict to the 
special case that 9' has no contribution besides these poles in order solve the recursion. We 
now have two parameters left in 9' . The distance of the poles from the real axis determines 
the exponential decay rate in the in the off-diagonal elements of the Hamiltonian and the 
strength of the residue 7 determines the pre-factor in front of the exponential. As is done 
in |Ha,Toj . we could get rid of the parameter by rescaling time, but we choose not to do 
so because tc plays a nontrivial role in optimal truncation and the error bounds obtained 
therein, and keeping this parameter will make things more transparent. 

We use (|43|) in order to determine the asymptotics of z„ . From Proposition |21 together 
with (|42() it is clear that ?;„ must go to zero dst ^ ±00. This fixes the constant of integration 
in H43|l ■ and we arrive at the linear two-step recursion 

zn+2{t) = (^z'^{t) + 9\t) J' 9'{s)zn{s) ds^ . (62) 

The fact that the recursion is linear will make its analysis simpler than the one of the nonlinear 
recursion in jHaJoj . We rewrite 9' as 

^'(t) = ?(/ + /) with fit)- 



tc t + itc 

For z„, we will make an Ansatz as a sum of powers of / and /. The reason for the success 
of this approach is the fact that this representation is stable under differentiation and inte- 
gration, and also under multiplication with 9' through the partial fraction expansion. More 
explicitly, we have 

Lemma 4. For each m> 1, 

m — 1 

6''Im(/™) = - 51 2-'=Im(/'"+i-'=), (63) 



/ m — 1 



6i'Re(/™) = - ( XI 2"'^Re(/"+i-'=) + 2-™6>' ) , (64) 



, /£ = 



Im(r)' = -!^Re(r+i), (65) 

777 

Re(/™)' = _lni(/'"+i). (66) 



— 2i ^ — 

Proof. We have f + f ^ = '^ff, and thus 



f'f = lf-\f + f) = lif' + f'-'f) 



and 

7 



ip-j = ^ I X f''+^-^^+^^ + 2""+V 1 ■ (67) 



Taking the complex conjugate of H67|) and adding it to resp. subtracting it from H67|) . we 
arrive at ^ and To prove (|65|l and H66() . it suffices to use that {f^)' = kf^^^ / (itc) 

along with the complex conjugate equation. □ 
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(n) 

Proposition 5. For each even n G N and j — 0, . . . ,n—l, let the numbers Uj be recursively 
defined through 

= 1, a^^^=0, (68) 



(n + 1); 



fc=0 m=0 



Jn+2) _ (n+2) (n+2) _ „ 

"ra — "n-1 J "n+1 ~ " • 



n-1 

^-^;r^E2"-''4"^Mr"^') (ne«en), (70) 



2/„ = 7^^^^) '2-M— > ar)Re(r-^) (n e«en) , (71) 



i=0 



n-l 



Proof. We proceed by induction. We have Xi = i6''/2 = i^Re(/), and thus by H38|l and 



Z2 = pT'i = -^Im(/2). 

This proves H7()|l for n = 2. Now suppose that 1)70(1 holds for some even n e N. Then by ((38|l 
and ((66(1 . 1(72(1 holds for n — 1. To prove 171(1 for the given n, we want to use ((39(1 . (163(1 and 
the induction hypothesis on z„ yield 

. , , n-l n-i-1 

*c j=o fc=0 



n-l 



J^E2-"ME«?MMr+^-'")- (73) 



-7 

Since ((73(1 only contains second or higher order powers of /, it is easy to integrate using l(66() . 
Let us write 



^^:^ = ^E«r- (74) 

Then by 1(66(1 we obtain 

yn = - f e'{s)zr.{s) ds = ^^^^ 2-"6(^) Re(/"'™) , 



proving ((71(1 for n. It remains to prove ((70(1 for n + 2. We want to use 1(6211 . and therefore 
we employ ((64(1 and our above calculations in order to get 

0'{t) f e'{s)zn{s)ds = 



15 



V^^^E^.I E 2-('=+^)Re(r+i-('=+^-)) + 2-"^?'U 

3=0 \ k=o ) 



(n-\)\ 



t 
By 



n+l 



'n-l / ^ \ \ /""I \ 



Now we sum the last two expressions, differentiate again and obtain 



(n-l) 

Zn+2 ^ -1 „+2 



,i=o \ fc=o / 



\k=0 



Comparing coefficients, this proves (|7UII for n + 2. 

(n) 

We now investigate the behavior of the coefficients as n — > oo. 
Proposition 6. Let a^""* be defined as in Proposition\^ 

(b) There exists Ci > such that for all n (z N 



n- 1 

(c) For each p > I there exists C2 > such that for all ri G N 

-5 1 (")l / ^2 

supp "'la; 'I < 



j>2 " n-l 

.(2) 



Froo/. (a) By a\;' ^1, and 

Comparing with the product representation of the sine function, 

^2 - 



(n+2) _ in) / ^ _ 2_ 



( 

3in(7ra;) = [] ~ ^) ' 

n— 1 ^ ' 



we arrive at (a). 

(b) Put a„ = (n - l)a^"'. Then by (jSHI), 
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thus for n — 1 > 7, we have 

|a«+2| < \an\ +7M ^ ~T ) ^axl^"^!, 

which shows (b). 

(c) Put Cj"^ — {n— \)p^^a^p\ and c^"^ = maxj>2 Icj""*!. We wiU show that the sequence c^") 
is bounded. We have 

\ fc=0 fc=l / / 

Now 

J k 



< r-^. (76) 



fe=2 m=2 

and 

We plug these results into H75|l and obtain 



p--V^< ^-^"^^^ < , \, . (77) 
•■^ — ' n — K n — J [n — J) Inp 



I („) / (n + 1 - j) (n + 1 - j) 7 V !_ 



\ n(n — 1) (n — — 1)2 n(n — 1) 

n (n — j)(p — 1)^ Inp 

By (a) and (b), Aq"^ and a^"'* are bounded. Taking the supremum over j > 2 above, we see 
that there exist constants Bi and B2 with 



n n{n — 1) / n 

hence 

^(«+2)_^(„) < i (^(^_2 + --^^c(") 

Now let n — 1 > _Bi . Then for c'"-* > B2 , the above inequality shows c*^""'"^-' < c*^"-' , while for 
c'"' < B2, c("+2) < c(") + B2/rt < ^2(1 + 1/n). Thus c(") is a bounded sequence. □ 

Remark 8. We will make no use of the fact that the logarithmic correction to the 1/rt-decay 
of the higher coefficients occurs only in the coefficient a^^\ We chose to include this in the 
statement of the preceding theorem anyway, because this gives some insight into the nature 
of the recursion and is not hard to prove. 

Remark 9. Numerical calculations of the first few thousand aj"-* suggest that (c) above 
continues to be true if we choose p = 1, but this seems to be much harder to prove. However, 
the estimate above is more than good enough for us. 
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Remark 10. The constants appearing in the proof of Proposition |51 (b) and (c) are not 
optimal, and could be improved by more careful arguments. This is unimportant for our 
purposes, and for the sake of brevity and readability we chose to use the simple estimates 
given. 

Corollary 3. Let 6^"^ be given by {74}) ■ Then for each p > I, there exists C3 > such that 

snpp ■'0' < 

J>0 



n — 1 

Proof. For j < n — 1, we have ri — 1 < j {n — j) , and thus Proposition IHl (c) gives 

,'3 / C2 p2 



in) 



< 



,(") I „(") 



n — 1 p — 1 



(p^-i - 1) < 



< P 



71-1 



< 



p — \ n — 1 n — 1 



□ 



Having good control over the coefficients a["^ , we can now derive relatively sharp estimates 



on the functions a;„, ?;„ and z„. Let us fix a < 1 and define 



Kit) 



(n - 1)" 



t + itc 



1 

V2 



n-2 



t + itc 



Obviously, for t < t^ the first function in the maximum above dominates, for t > tc the 
second one docs. 

For families of functions gn{t),Gn{t) we write 

if there exists C > such that \gn{t)\ < C\Gn{t)\ for all n e N and aU t e R. 
Theorem 3. For n> 1 and a <1, we have 

.{n-l)\ ( 2sin(77r/2) 



Xn{t) = 

yn{t) - 

Zn{t) = 



Re 1 - i- 



o{K{t)) 



(n-l)! 



o{K{t)), 



{n- 1)1 / 2sin(77r/2) 



Im I I 1 - i— 



(78) 
(79) 
(80) 



Proof. With the definition of / and Proposition IHl (a) we get 



4")lm(/") 



2sin(77r/2)^ 
7r7 



Im 1 - i- 



t + itc 



when n is even, and a similar formula for aQ""''^^Re(/") when n is odd. This covers the j = 
terms in (|72f) and jTUJ). For the remaining terms, let 



An) 

'"7 

-,(«) 



if n is even, 



na^- ' /{n — j) if n is odd. 
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Now n/{n — j) < j for j < n, and thus by Proposition (b) and (c) for each p > 1 we can 
find C > such that 

(") / • 7 ^ 

c) <]p'- 



{n - 1)" 

for all j > 1. (For j > 2, we may even choose a — 1, but we will not exploit this.) For 
\t\ < tc, we have \tc/{t + iQl'^ < 2-''/2, so we get 



E 

J=2 



1 ^(«) 



21 \t + it. 



< 



c 



in-1) 



i=2 



V2j 



) 


tc 




t + itc 



If we choose p < \/2, the sum on the right hand sided is bounded uniformly in n. Combining 
this with our above calculations, H78|l and 1)80(1 are proved for |t| < tc. For \t\ > tc, we have 
\tc/{t + itc)\ < 1/V2, and thus 



n-2 (n) 



E 

J=2 



2J \t + itr 



< 



c 



{n - ly 



t + itc 



2 n-2 



^^■(f) 172 



J=2 



«-2-j 



If we choose again p < \/2, the sum on the right hand side is bounded by (7(1/^2)" ^ 
uniformly in n. For the term with j = n — 1, this does not work since then n — 2 — j < 0. 

(n) 

But for n even, this term vanishes since then c^li = 0, and for n odd, it equals 



„(") 



2" 



-Re 



itr 



t + itr 



(n+1) 



2" t2 + t^r 



tl ^ c 



- 1 V\/2 



^ ri-2 






< + itc 



This proves lf75|l and l|5n|) for |i| > tc- The proof of O is similar and uses Corollary|3| □ 



5 Optimal truncation 

By the results of the previous section tt^ grows like (fc — l)\/tc- Hence, the sum tt^"^ = 
^^^g e'^TTfc does not converge to an exactly equivariant projection tt*-""^ as n — > oo. This is 
the reason why we see exponentially small transitions. The basis in which these transitions 
develop smoothly is the optimal superadiabatic basis: since we cannot go all the way to 
infinity with n, we fix e and choose n — n(e) such that the off-diagonal elements in jSJ 
become minimal. Using Stirling's formula and 1(78(1 resp. I(8()() . it is easy to see that the place 
to truncate is at n[e) — tc/s. This is in general not a natural number, but we will find 
that a change of n which is of order one does not change the results. Before we go into more 
details, we need a preliminary result. 

Lemma 5. Uniformly in x £ [0, 1] and for k > 0, we have 

(l + x)-'' =.e-'=^+e^'=^/^C'(i). 
Proof. We start with the equality 

(1 + 2;)-'= - e"*^^ = e"*^^ |^efc(:i;-in(i+x)) _ ^'^ _ ^g^-j 

At first consider x > \f\Jk. There we use the inequality (a; — ln(l + x)) < a;/3, valid for 
< a; < 1, in (|FT|l and obtain 

1(1 + a;)-*^ - e"*^^! < e-'^'^ (e'='=/^ " ^) " '"'^^^ {^-^''1^ - e-'^^^/^j . 
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For X > \/l/k, the term in the last bracket above is 0{l/k), and we are done in this case. 
For X < -\/l/fc, we use (x — ln(l + x)) < and rearrange H81|l to get 

To find out where f{x,k) is maximal, we calculate 

±fix,k) = ^e~'^/' {l + e'^'/\2x-l)). 

The derivative is zero exactly at the solutions of the equation 

ln(l - 2a;)/.T2 = -A:/2. (82) 

Now ln(l — 2x)/a;^ = —2/x + R{x), where R{x) is a power series in a;, convergent for a; < 1/2. 
Thus for x < y/l/k and k sufficiently large, there exists exactly one solution x*{k) of l|82(l . 
and x*{k) < C/k uniformly in k for some C > 0. Since 3^/(2;, fc) > for a; < f{x, k) 

has a maximum at x* (k) . Thus 

fix, k) < f{x*{k),k) < e-^/2fc _ 1 ^ 0(1/^) 
for X < \f\Jk, and the claim is proved. □ 
Lemma El immediately yields 

(l + ^)"'=e-(l + 0(i)) (83) 

uniformly on compact intervals of a by taking x = ajk. 

We now turn to the proof of Theorem ^ which we deduce from Theorems El and O As 
stated already in we will use 

ne = --l + cre, (84) 

e 

where G [0, 2[ is such that is an even integer. The advantage of this convention 
about (Te is that now the off-diagonal components in (|60|) are always given by £"=+^x„^+i 
since z„+i — for even n. Of course we could as well consider the asymptotic behavior of 
e^+^Zn+i for odd n and one would expect to end up with the same result. However, it is 
obvious from H78|l and H8U|) that Xn+\ is purely imaginary and Zn+i is real at leading order. 
Thus the large n asymptotics of the off-diagonal component of the effective Hamiltonian do 
depend on whether we consider even or odd superadiabatic bases. On the other hand, the 
asymptotics of the propagator must be independent of the exact choice of basis. We will 
discuss this point after giving the proof of Theorem phased on the above convention. 

Proof of Theorem^ We want to apply Theorem |2 and thus have to check that ^,rj and ( 
defined in (|54|I - H56|I together with their derivatives are 0{ed'). From Proposition [3 together 
with Proposition El we infer that there exists C > such that |a;fc(i)| < C9'{t){k — 1)!/^^ for 
each k. The same is true for y„ and z„. Using the differential equations we find 

that there is C" > with K(i)| < C'9' {t)n\/t]i+'^ . This means that 

n 

W{t)\ < eC'e'{t) ^ eh-^-^kl e'^-i, 
fc=i 
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with similar expressions for the other quantities. Now taking e = tcjijie ~ <^e\ we find 

(fc + 1)! 2 3! 



fe=0 



fc=0 



rig — (Tj (rie — cr^ 



Each of the + 1 terms in the sum above is bounded by const/ (n^ — a^) except the first 
which is 1. This shows 

and Theorem|21gives (O with c"-(t) = e"-+^a;„^+i(t)(l + 0(e6''(t)). Recall that z„^+i(t) = 
due to our convention. It remains to determine the leading order asymptotics of e"=+^a;„^+i. 
For convenience of the reader let us rewrite (|78|l as 



£"= + ^X„. + i(t) =i 



4-rie + l 



2sin(77r/2) 



Re 



(85) 



Lemma 6. With \<iJ^ , we have 



Cc V 

Proof. Stirling's formula for (n + 1)1 implies 

n+l 



1 /n+l 



n + 1 \ c 
Together with this yields 



VnTTv2^ 1 + o 



He + 1 



2tt 



n, + 



C^+^e" 



27r£ 



Finally, 



tc + eae 



2ne / eae 
IH 



-1/2 



27r£ 



□ 



Lemma El takes care of the first factor in (|85|l . Turning to the terms inside the square 
brackets in H85|l . let us first note that for \t\ > tc, both terms are 

0(2-(»^"i)/2/(i + t2)) = 0(exp(-ie ln2/(2e))/(l + t^)) , 

proving the theorem in this case. For \t\ < tc, we investigate the modulus and the phase 
separately. Let < (3 < 1. From Lemma |5l it follows that 



1 + 



"e + l 



= 1 
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For \t\ > e'^/^, exp(-t^/(2tc£)) = 0{e cxp{-P /{itcs))). Thus neither the prefactor hivolving 
CTg above nor the phase play any role in this region. For \t\ < e^^"^, = l+0{a^e^) 

and therefore 



1 + 



The same reasoning applies to R^^^i and gives 



Turning to the phase in the region |t| < e^^'^, we find 



e 



exp 
exp 



(i (l - 3^ + ^)) (l + Ois""/'-') + Oie'^^')) . 



Now we just have to collect all the pieces and add the complex conjugate. □ 

Let us now see what of the above would have changed for odd. Then Xn^+i = 0, and 
l|Sn|) together with Lemma |S1 and yields 

c^-W = -e"^+iz„,+i(i)(l + 0(e(?')) (86) 



At first, this looks like an important difference, since now the off-diagonal elements in the 
transformed Hamiltonian are purely real- valued in leading order, while in the other case they 
were purely imaginary. However, in the computation of the propagator, another factor of 
exp(±ii/e) from the dynamical phase appears, cf. H87|l . At leading order only the resonant 
term of the Hamiltonian survives, which is the same for odd and even rig. 

6 First order perturbation in the optimal superadiabatic 
basis 

In this section we prove Corollary ^ Since we use standard first order perturbation theory, 
we stay sketchy in some parts. After splitting H"'{t), see as 

4 ) + Veit) H, + V,{t) , 

Dyson expansion in the interaction picture (cf. |R,eSij . Thm. X.69) yields 

K^'it,s) = e-^ fid- - / e^ I4(r)e — ^drje-r' 

\ 0(ee-^) 
Thus we only need to evaluate the integral 



A(t,s) 



-- / e e V^frje ^ dr = -- / ^ , , 



0(e) 
0{e) 



Ait,s). 
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Inserting H14|l and using IjlSfl gives 
J e^c"=(T)dT = 



sin — ^ e"~ / e~ e"~ e " ^ + e " ^ Mr 



+ o(£"e-^A(i,s)j (87) 

for each a < 1. Now we replace the exponentials e ^^*<: by 1 ± — Using 

je"'' — 1 — iLp\ < ip'^, we conclude that the resulting error is bounded by a constant times 

£ 2e'- J e 2=tc ( _ + _ + r^^j dr = 0(ee"~) . 

Hence we obtain 



W = V^sin(^)e-^y^e-^(l + |J-l^+e^(l-i^ + l^)) 



with a < 1, where the first summand in the integrand gives rise to the explicit term in p7|) . 
The remaining terms can be integrated explicitly as well, most conveniently using Maple or 
Mathematica. They are all of order 0{y/ee~i' A{t, s)) uniformly in t and s resp. of order 
0(£"e~^ A(t, s)) for \t\ and \s\ larger than e'^ for some P < \- To illustrate the reasoning 
note that 

e 2etc r dr = rf^ ( e ^etc - e ^et^ j . 

This is uniformly of order ©(e), but of order 0{e~^^'^ ^ ) for \t\ and \s\ larger than . Finally 
we emphasize that we could get the next order corrections to (I17|l explicitly by evaluating 
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